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We consider the Neumann Sturm-Liouville problem denned on trees such that 
the ratios of lengths of edges are not necessarily rational. It is shown that the 



Abstract 



potential function of the Sturm-Liouville operator must be zero if the spectrum is 
equal to that for zero potential. This extends previous results and gives an Am- 
barzumyan theorem for the Neumann Sturm-Liouville problem on trees. To prove 
this, we compute approximated eigenvalues for zero potential by using a general- 
ized pigeon hole argument, and make use of recursive formulas for characteristic 
functions. 
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1 Introduction 



In this paper we study an eigenvalue problem for a Neumann Sturm-Liouville op- 
erator defined on a metric tree (a connected graph with no cycles) T = {V,E}, where 
V = {vj : j = 0, . . . , J} is the set of all vertices and E = {7^ : i = 1, ...,/} is the set 
of all edges with lengths a* G (0, 00). We let dT be the set of all pendant (boundary) 
vertices. Any edge connected to a pendant vertex is called a boundary edge. For any 
internal vertex v, we let J(v) be the set of all indices % such that the edge 7, is connected 
to v. The degree of v is defined as |/(v)|. We also choose an arbitrary internal vertex 
v to be the root. Then we assign local coordinates to the edges such that for any 
edge oti G E, the endpoint of 7$ further from the root has local coordinates 0, while 
the endpoint closer to the root has local coordinate a%. Thus all pendant vertices have 
local coordinate 0, while the root vo has local coordinate dj corresponding to the edge 
7i connected to it. Thus the Neumann Sturm-Liouville problem can be expressed as 
functions y^s defined on the edges Tj's satisfying 

- y" + qi(x)yi = Xyi 0<x<ai, 1 = 1,2,...,/, (1.1) 

where each g« is the potential function defined on % (i = 1, . . . , I), and, 

(Al) y'i(0) = whenever 7^ is a boundary edge. 

Also continuity and Kirchhoff conditions are imposed at each internal vertex. At 
any internal vertex v other than the root, there are some incoming edges 73 's and one 
outgoing edge 7^. The local coordinates of v is with respect to 7^ but with respect 
to the other 7j's. 

(A2) The continuity and Kirchhoff conditions at v are defined as 

l/*(0) = yM, and ^(0) = ^^(a^), 

i 

wherever % is such that % G /(v) but i 7^ k. 
At the root vo, all the connecting edges are incoming. Hence 
(A3) the continuity and Kirchhoff conditions are defined as: for any i, k G 

Vk(ak) = i/i(oi), and ^ Vii^) = °» 

ie/(v ) 

The above formulation of the Sturm-Liouville problem on F is essentially the same 
as in [TU]. Note that the Neumann boundary conditions (Al) can be viewed as a special 
case of (A2), for the degree of a pendant vertex is 1, hence continuity condition is empty 
and right hand side of the Kirchhoff condition vanishes. It is well-known that the above 
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problem has a discrete spectrum. A real number A is an eigenvalue of the above problem 
if it has a nontrivial solution (y\, . . . , yj), i.e., at least one of y^s is nontrivial. We write 
Q — (qx, . . . , qj), and let a(Q) be the set of eigenvalues for the vector potential function 
Q. In particular, cr(0) denotes the set of eigenvalues for Q = 0, i.e., = for all 
i = l,...,J. 

In the simplest case 1 = 1, the above problem is reduced to a usual Sturm-Liouville 
problem on a finite interval (0,ai) with the Neumann boundary condition. In this case, 
in 1929, Ambarzumyan pQ showed that a(Q) = er(0) implies Q = almost everywhere. 
This seems to be the first example of an inverse problem where the potential function 
can be determined uniquely by the spectrum, without any additional data. Later, Chern 
and Shen [I] extended the result to vectorial Sturm-Liouville systems. In the case of 
the Dirichlet boundary condition, it was shown by Chern et al. [3] that the potential 
function must be identically equal to if an additional condition 

qi(x) cos( — x) dx = 
ai 

is imposed. Recently, the Ambarzumyan problem for periodic boundary conditions was 
studied by Yang et al. [14J. They showed that for a vectorial Sturm-Liouville system of 
dimension d with the periodic boundary condition, if the eigenvalues are (2mr) 2 with 
multiplicities 2d, then the vector potential must be 0. Thus the Ambarzumyan theorem, 
originally specified for the Neumann boundary condition, can be generalized to several 
Ambarzumyan problems with different boundary conditions. 

The aim of this paper is to study the Ambarzumyan problem for the Neumann 
Sturm-Liouville problems defined on trees. In this direction, Pivovarchik [11] showed 
that when F consists of three edges with equal length and one triple junction, an analogue 



of Ambarzumyan theorem is valid for the problem ( 1.1 ) with (Al)~(A3). Later, Carlson 
and Pivovarchik [2] extended the result to any trees such that {cii/L} are all rational 
numbers, where L is the total length of T given by 



l 

a. 

i=l 

In this case, we can find infinitely many eigenvalues explicitly given by 

(^p7r) 2 Ga(0), m= 0,1,2,... (1.2) 

for some integer mo, which makes the analysis much easier than a more general case. As 
for the Dirichlet problem on a star-shaped graph, we refer to a recent paper by Hung et 
al [5] . See also [HI [TJ EJ |9j QUI El [13] for related results on the Sturm-Liouville problems 
on graphs. 

In this paper we consider the Sturm-Liouville problem on trees such that {a^/ L} are 
not necessarily rational. The following theorem is a main result of this paper, which 
gives a solution to the Ambarzumyan problem on general trees. 
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Theorem 1.1. For the Neumann Sturm-Liouville operator defined on T, a(Q) 
implies Q = almost everywhere. 



<t(0) 



When we deal with general trees, we encounter two kinds of difficulty. The first one 
is that we must handle trees with arbitrary number of edges. Moreover, even if the 
number of edges is given, there are various trees with different topology. In order to 
handle all trees, we shall derive a recursive formula for characteristic functions whose 
zeros are the square roots of eigenvalues. The second one is that we may not have explicit 



eigenvalues as (1.2). To overcome the difficulty, we approximate {ai/L} precisely by 
rational numbers at the same time by applying a (generalized) pigeon hole argument. 

In Section 2, we shall study direct problems and derive recursive formulas for char- 
acteristic functions. In Section 3, we compute the expansion of characteristic functions. 
In Section 4, we present a key lemma for the approximation of eigenvalues. Finally 
Section 5 is devoted to the proof of Theorem 



1.1 



Hereafter in this paper, we shall let Vj be a pendant vertex which is an endpoint of 
the edge 71, while v 2 is another vertex at the other end. And without loss of generality, 
we assume that 7 2 is connected to v 2 , and is closer to v than 71. 



2 Recursive formula for characteristic functions 

Let p > 0, and let y = Ci(x; p) and y = Si(x; p) be (linearly independent) solutions of 

- y" + qi(x)y = p 2 y, < x < a { (2.1) 

with the initial conditions 



y(o) 



and 



y'(0) = 0, 
1. 



1/(0) = 0, y'(0) 

respectively. We sometimes write Cj(x; p) as Cj(x) and Si(x; p) as Si(x) just for sim- 
plicity. For each edge 7$, we write y,i as 

Hi = Aid(x) + BiSi(x). 

Then from (Al)~(A3), we have a system of linear equations for the unknowns {A} and 
{Bi}. Thus we may express the coefficient matrix of the system of linear equations as 



$n(p) 



C((0) S[(0) •■• 

d(ai) Si(aa) -C 2 (0) -S 2 (0) ••• 

Ci(ai) S[(ai) -C 2 (0) -S' 2 (0) * ■■■ * 

* * * • • • * 



(2.2) 
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where the first row corresponds to the Neumann boundary condition at Vi G T, and 
the second and third rows correspond to the continuity condition and the Kirchhoff 
condition at V2. Note that in the first row, we include the term S[(0) although the 
corresponding coefficient B\ = 0. We do this in order to have a systematic form of $at 
during reductions, as we shall see later. We define a characteristic function by 

<p N (p) := det$ N (p), 

so that A = p 2 G u{Q) if and only if <Pn{p) — 0. We note that the characteristic function 
depends on the orientation of edges and how to express the coefficient matrix, but the 
set of zeros of ^Pn{p) does not depend on them. 

Next we introduce another eigenvalue problem by replacing (Al) with the following 
condition: 



(A4) If V! is an endpoint of any boundary edge 7^ the solution of (1.1 ) satisfies the zero 
Dirichlet boundary condition 



2/i(0) = (or yfa)) = 

At other boundary vertices, the solution satisfies the homogeneous Neumann 
boundary condition as in (Al). 



Hereafter, we call (1.1) with (A2)~(A4) the Dirichlet-Neumann problem. If we impose 
the zero Dirichlet condition at 01(71) G dT, then Ax = 0. Thus we may express the 
corresponding coefficient matrix by 



®d{p) 



Ci(0) S'i(O) 

Cx(ai) S 1 (a 1 ) -C 2 (0) -S 2 (0) 

C[{a x ) S[( ai ) -C' 2 (0) -S' 2 (0) 

* * 











* 
* 



Then we define a characteristic function for the Dirichlet-Neumann problem by 

(f D (p) := det$£>(p). 

Again, the set of zeros of (Pd{p) does not depend on the orientation of edges and how 
to express the coefficient matrix. Our interest will be only in zeros of the characteristic 
functions, and hence the non-uniqueness of characteristic functions will not affect the 
following argument. 

For a general tree, it is not easy to express explicitly the characteristic functions (p^ 
and Instead, we may compute these functions recursively as follows. Let T be a 
subtree of Y obtained by removing 71. We denote by (f>N the corresponding characteristic 
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fig.l 



function of the problem with Neumann condition at v 2 in case J(v) = {1, 2}. However if 
the degree of v 2 is greater than 2, then we take the continuity and Kirchhoff conditions 
at v 2 instead. Thus we say (pN the characteristic function for a Neumann/Kirchhoff 

problem on V. (see fig. 1) 

If we remove 71 and replace the matching conditions (A2) and (A3) by yj = at the 
vertex v 2 for any j € (v 2 ), we have subtrees of T on which the Dirichlet- Neumann prob- 
lems are defined. For each subtree, a characteristic function of the Dirichlet-Neumann 
problem is defined as above. We denote by <Pd(p) the product of these characteristic 
functions. 

For trees with two or more edges, we have the following recursive formulas. 

Lemma 2.1. Assume that I > 2. Then the characteristic functions have the following 
properties: 

(a) (p N (p) = Ci(a 1 )(p N (p) - C[(a 1 )(p D (p). 

(b) (p D (p) = -S , i(ai)(^ A r(p) + S' 1 (a 1 )(p D (p). 

Proof. Since C[(0) = and 5^(0) = 1, by expansions with respect to the first row and 
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the first column, we have 



det $at(p) = - det 



C[( ai ) 






-C 2 (0) -S 2 (0) 
-C 2 (0) -S' 2 (0) * 





* 



= Ci(ai) det $jv(p) - Ci(ai) det 
where l>7v and $d are (21 — 2) x (2/ — 2) matrices given by 



$jv(p) 



C 2 (0) S 2 (0) * 



* * 



and 



C 2 (0) 5 2 (0) 





* 



* * 



respectively. Noting that &n describes conditions on the Neumann/Kirchhoff problem 
for T with 7x removed, we have 

det$Ar(p) = (p N {p)- 

Similarly, since $£> describes conditions on the Dirichlet-Neumann problems for subtrees 
of T obtained by removing 71, we have 

det$Ar(p) = (f> D (p)- 

Thus the proof of (a) is completed. 

Next, let us consider the Dirichlet-Neumann problem. In this case, we have the 
expansion, since Ci(0) = 1 and Si(0) = 0, 



det $d(p) = det 



Si(ai) -C 2 (0) -S 2 (0) ••• 
S[(ai) -C 2 (0) -S' 2 (0) * ••• * 
* * * • • • * 



* * * • • 

-51(01) det $ N (p) + S[( ai ) det $ D (p). 



This proves (b). 



□ 
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Remark 2.2. More general recursive formulas were obtained in a recent paper by Law 
and Pivovarchik [W{ . Interested readers might like to read a spectral determinant ap- 
proach to the same formulas [Wj . 

Now we compute the characteristic function of a general tree as follows. Given a tree 
with two or more edges, we remove one of the edges of F and use the recursive formulas. 
Repeating this procedure, we will reach to problems on single edges. For a single edge 
7i, we may define its characteristic functions by 

<Pn(j>) ■= -C'i(ai), <Pd ■■= S-idi). (2.3) 

Thus we can express <Pn(p) an d <Pd(p) as polynomials of {Cj(aj)} and {^(cij)}. 

Next, we consider the zero potential Q = 0, and denote by ip n (p) and iPd{p) the corre- 
sponding characteristic functions of the Neumann/Kirchhoff problem and the Dirichlet- 
Neumann problem, respectively. Similarly, we denote by ipN and ^ be characteristic 
functions for the Neumann and Dirichlet-Neumann problems for T with ji removed. 



For a tree with a single edge 7j, to be consistent with (2.3), we define its characteristic 
functions by 

i) N {p) ■= sm(pai), $d{p) ■= cos(paj). (2.4) 
For trees with two or more edges, the characteristic functions can be computed by using 



the following recursive formulas repeatedly and (2.4) 



Lemma 2.3. Assume I > 2. Then the characteristic functions for Q = have the 
following properties: 

(a) ip N (p) = cos(pai)V>iv(p) + sin (pa x)$ D {p)- 

(b) ip D (p) = -sin(pai)^Ar(p) + cos(pai)^ D (p). 



Proof For qi(x) = 0, the solutions of ( |2.1 ) are given by Ci(x; p) = cos(px) and Si(x; p) 



sin(px) . Then the above recursive formulas can be obtained in the same way as Lemma 2. 1 



□ 



3 Expansion of characteristic functions 

In this section we show the following asymptotic formulas for Cj(aj; p) and Si(ai \ p) as 
p — > oo. 

Lemma 3.1. As p — > oo, one has 

Ci(au p) = cos(pai) + p~ l Ki sin(paj) + o(p _1 ), 
C-(a;;p) = -psin(pai) + K, t cos(pai) + o(l), 
Si(ai] p) = p' 1 sin(paj) - p~ 2 Ki cos(paj) + o(/T 2 ), 
S[(ai] p) = cos(pai) + p~ x Ki sin(pai) + o(p _1 ), 



S 



where 



Hence 
and 



1 f ai 
Ki := 2 J qi ^ dx - 

Proof. By the variation of constants method, it is easy to show that Ci(x) satisfies 

1 f x 

CAx) = cos(prr) H — / sm(p(x — tUqAtjCAt) dt. 
P Jo 

Ci(x) = cos(px) + 0(p- r ) 

Ci(x) = cos(px) H — / sin(p(x — t))qi(t) cos(pt) dt + 0(p~ 2 ) 
P Jo 
1 f x 

= cos(px) + — / (sm(px) + sin(p(:r - 2t))qi(t) dt + 0(p~ 2 ) 
^P Jo 

= cos(px) + sm ^ I q^dt + oip- 1 ). 
^P Jo 

Furthermore, 

C[(x) = -psin(px) + / cos(p(x - t))qi(t)Ci(t) dt 

Jo 

= —p sin(pa;) + / cos(p(x — t)) cos(pt)qi(t) dt + 0(p^ 1 ) 
Jo 

1 f x 

= — psm(px) + - / ( cos(prr) + cos(p(:r — 2t)))gj(t) dt + 0(p _1 ) 

4 .In 



COS( OX) f x 

-psin(px) + — f-t- / qi (t)dt + o{l). 
^ Jo 



Evaluating at x = a j, we obtain the asymptotic formulas for Cj(aj) and C^(aj). 
The asymptotic formulas for S^Oj) and S^(aj) can be derived similarly by using 

= ^m(px) + ir _ t ^ q ,^ s ^ dt 

P P Jo 

We omit the details. □ 

Lemma 3.2. The characteristic functions ifN and <p D have the following properties: 

i 

(a) (fi N (p) = P^n(p) ~ (^2 K i)^D{p) + o(l) as p ->• oo. 

i=i 

/ 

(b) (po(p) =i/>d(p) + P~ l (^2K^i) N (p) + o(p~ 1 ) asp-; 

i=i 



Proof. If T consists of a single edge, then by Lemma 3.1, we have 

¥n{p) = -Ci( a i) = psin(pai) - K 1 cos(pa 1 ) + o(l), 
<Pd(p) = Si (ax) = cos(pai) + p~ 1 K 1 sin(pai) + o(p^] 



so that 



<Pn(p) = p4>n(p) - Ki^ D (p) + o(l), 
<Pd(p) = iPd(p) + p' 1 K 1 i) N (p) + o(p^). 



Hence (a) and (b) hold in this case. 



Suppose now that (a) and (b) hold for any subtree of V. Then by Lemmas 2.1 and 



2.3, we have 



Vn{p) = Ci{a l )ip N (p) - C[{a 1 )(p D {p) 

i 

= {cos(pa!) +p- 1 iT 1 sin(pa 1 ) + o(l)} j/«Mp) - (j^ K^jip D (p) + o(l)\ 

i=2 

I 

- | -psin(pai) + K x cos(pai) + o(l)| {^(p) + ( ^ K^jip N {p) + o(p~ l 

i=2 

= pcos(pa 1 )^ Ar (p) +iT 1 sin(pai)'0 i v(p) - cos(pai) [ ^ K^)ip D (p) 

i=2 

I 

+ ps\n{pa x )$ D {p) -K x cos(pai)^ D (p) + sin(pai) ( ^ Kj)ij) N {p) + o(l) 
= p^(p) - #iV>d(p) - ( K i)Mp) +o(l), 



i=2 



i=2 



1 

p^(p)-(E^)^ d(p) + o(1) - 

i=l 
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Similarly, 

<Pd(p) = -5 , i(ai)(^ A r(p) + S' 1 (a 1 )(p D (p) 

= -{p 1 sin(pai) - p~ 2 K 1 cos(pai) + o(p~ 2 ) j jp^Ar(p) - ( ^^^(p) + o(l)| 

+ jcos(pai) +p~ 1 K 1 sin(pai) + o(l)} {^d(p) + P _1 ( ^ K^ N {p) + o(p- 1 )} 

i=2 

= - sin(pa 1 )^ Ar (p) +P _1 j #1 cos(pa 1 )i) N (p) + sin(pa 1 ) f y] K^ D (p) \ + o(p _1 ) 

^ i=2 J 

+ cos(pai)^ D (p) + p _1 jfTi sin(pai)^ D (p) + cos(pai) f ^ ifr) ^jy(p) | + o(l) 
= Vd(p) +p- 1 {^i^(p) + (^^)^(p)} +o(p- 1 ) 

^ i=2 J 

= Mp) + p - i (J2 k ^^ + 



i=2 



i=l 



Hence the assertion holds for the whole tree. Thus by induction, the proof is complete. 

□ 



4 Approximation of eigenvalues 

In this section we compute approximate eigenvalues for T with Q — 0. We begin 
with the following lemma. 

Lemma 4.1. There exist infinite sequences of natural numbers {m n } and {k itn } {i = 
!,...,/) such that m n — > oo as n -> oo and 



< m," 



-i-i// 



for all i — 1,2, I and n = 1,2,.... 

Proof. If a«/L = 1, . . . ,1) are all rational numbers, we can find natural numbers p 
and qi such that 

a i Qi , „ 

y = -, i = l,2,. ..,1. 
L p 



Then we may take 



h,n = nqi, m n = np. 
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Assume that not all of a>i/L (z = 1, . . . , I) are rational. We consider the /-dimensional 
unit cube 

{x = (x±, X2, ■ ■ ■ i xi) G R 1 : < Xi < 1 for i — 1,2, ... , I}, 

and divide it into n 1 smaller cubes with edge length 1/n. We shall approximate the 
vector (cii/L, . . . ,cii/L) by a rational point in the unit cube as follows. For each p = 
0,1,2,..., n 1 , there exists g ; G N such that 

o< P j- qi <i. 

Then by the pigeonhole principle, there exist vectors U\ and U 2 of the form 

a 7 



that fall into the same small cube. Then for m n = p 2 — pi < n 1 and h^ n = — q%,i, we 
have 



< 



n 



so that 



(Xi k>i,i 

L m r 



1 1 

< < 



V m 



1+1/r 



The proof is complete since from (4.1), m n — > 00 as n — > 00. 



(4.1) 



□ 



i=l 



1=1 



^2 ^i,n 



< 



n 



Remark 4.2. 

(a) Since 

the equality 
holds for all n > I. 

(b) Using X][=i( a «/-^) = 1> we ma V apply the pigeonhole principle on the (I — 1)- 

dimensional hyperplane 

{x = (xi, x 2 , ■ ■ ■ , Xi) E K 1 : xi + x 2 H h xj = 1}. 



5> 



i=l 



(4.2) 



TTjen we can improve the result in Lemma 4J_ to 

a,i hi 



'i,n 

L m n 
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Let {m n } be the sequence given in Lemma 4.1 We define a sequence {p„} 



2m n 7T 



n = l,2,..., 



and compute the values of ip^, ipr>, and their derivatives at p = p n as follows. 
Lemma 4.3. The characteristic functions ipN and ipo have the following properties: 

(a) i> N (Hn) = 0(pn l/I ) and — ^-(p n ) = L + 0(p~ l/I ) as n ->■ oo. 

dp 

(b) i) D {^n) = 1 + 0(pn 1/I ) and — ^(p n ) = Ofp" 1 ^) as n ->■ oo. 

dp 



Proof. By Lemma 2.3 and 



Pn^i = — ?Oi = 2fcj,„vr + 0(p n 1/7 ), 



we have 

ipN(Hn) = cos(p n ai)i) N {p n ) + sin(/i n ai)^£)(/i n ) = $ N {n n ) + 0(fj^ 1/I )^ D (n n ) > 

and 

= -sin(/i n Oi)^2vW + cos(/i n ai)^ D (/i n ) = ^o(/i n ) + O^" 1 ^)^/^)- 
Using these equalities repeatedly on ( |2.4 ), we obtain 

V^On) = 0(/i^ 1/7 ), ^7v(/in) = 0(^ 1/7 ), 

and 

^ D (/i n ) = i + o(p; 1/7 ), ?Mp n ) = i + o^ 1 / 1 ). 

Also, we have 



dp 



d 

dp 



cos(pai)^Ar + sin(pai)-0£) Hp 



d^ 



-ai sin(p n ai)^v(Pn) + cos(p n ai)^^(p n ) 



dp 

d?/> 



+ aiCos(p„ai)^D(Pn) + sin(p n ai)^p^(p n ) 

dp 



N 



dip. 
dp 



(Pn) + a 1 + 0(p, 
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and 



dtp 



D 



dp 



(fa 



d 
dp 



- sin(pai)^jv + cos(pai)ip D \(n n 



-ax cos(p n ai)%l) N (p n ) ~ sin(/x n ai 



dip 



N 



dp 

dp 



dip 



dp 



ai s\n(ii n ai)ip D (ii n ) + cos(/j, n ai)— r —(jj, n ) 

D (^n)+0(^ 1 ). 



Using these equalities repeatedly and (2.4), we obtain 

1 



dip 



N 



dp 



i=l 



and 



The proof is complete. 



dip 



D 



dp 



W = o( A ^ 1 /'). 



□ 



As an immediate consequence of Lemma 4.3, we have the following result concerning 
the location of zeros of ipN- 

Lemma 4.4. There exists a sequence of positive numbers {p n } such that ipN{pn) — 
for n = 1, 2, . . . and p n = p„ + 0(jJLn 1 ^) as n — >■ oo. 



Proof. By Lemma 4.3, there exists M > in dep endent of n such that \ipN(p n )\ < 
M/in. 1 ^ 1 for all large n. On the other hand, by (2.4) and Lemma 2.3, iPn{p) and iPd(p) 
must be some polynomials of sin(paj) and cos(paj) (i = 1,2,...,/) with constant co- 
efficients. This implies that \{d 2 / dp 2 )ip^\ is uniformly bounded in p > 0. Hence by 
Lemma 4.3, there exists 5 > independent of n such that 

—ipNip) > 77 for p G (p n - 5, /!„ + 5). 
dp 2 

Therefore iPn(p) must vanish at some p n such that 

|Pn-Pn| < {2M/L)^7 n 1/I . 

This completes the proof. □ 
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5 Proof of Theorem 



1.1 



In this section we complete the proof of Theorem LI By Lemma tp^ satisfies 

i 

<Pn(p) = Mv(p) ~ ( E ^i)V'o(p) + o(l) as p ^ oo. 

Here, setting p = p n , we have y?7v(p n ) = = ipN(Pn) by = cr(0). Also, by 

Lemmas 4.3 (b) and 4.4, we have 

MPn) = l + 0(p- 1 / 1 ) 



as n — )■ oo. 



Thus we obtain 



(E^){ 1+o(p " 1/7) } +o(1) = - 

i=i 



Letting n — >■ oo, we conclude 



i=l 



i=l 



(5.1) 



On the other hand, since A = is the first eigenvalue, it follows from the variational 
principle that the inequality 



i=l 



> 



1=1 



E / (Vi) 2 dx 



holds for any test function in iJ 1 (r). By (5.1), the infimum is attained by Ui{x) = 1 
for i = 1, . . . , I. This implies that the constant vector (1, . . . , 1) is an eigenfunction 
associated with A = 0. Then by simple substitution, we conclude that qi = a.e. on 
[0, di], i — 1, 2, . . . , /. Thus the proof is complete. □ 
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